Depletion potential in the infinite dilution limit 
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The depletion force and depletion potential between two in principle unequal "big" hard spheres 
embedded in a multicomponent mixture of "small" hard spheres are computed using the Rational 
Function Approximation method for the structural properties of hard-sphere mixtures [S. B. Yuste, 
A. Santos, and M. Lopez de Haro, J. Chem. Phys. 108, 3683 (1998)]. The cases of equal solute 
particles and of one big particle and a hard planar wall in a background monodisperse hard-sphere 
fluid are explicitly analyzed. An improvement over the performance of the Percus-Yevick theory 
and good agreement with available simulation results are found. 



I. INTRODUCTION 

Excluded volume interactions in hard-sphere mixtures 
are interesting for a number of reasons. For one thing, 
colloidal systems are often modeled as mixtures of dis- 
similar hard spheres and, experimentally, it is known that 
this interaction plays an important role in the observed 
behavior. Moreover, due only to the existence of a dis- 
parate size ratio between solute and solvent, the (purely 
entropic) effect on the solvent mediated interaction forces 
between solute particles in a hard-sphere suspension may 
be quite dramatic. Take for instance the case of two (not 
necessarily equal) big hard spheres immersed in a fluid of 
small hard spheres. When the distance between the two 
big spheres is less than the diameter of the small spheres, 
the latter may not got into the gap. This depletion effect 
induces an imbalance in the local pressure leading ulti- 
mately to an effective attraction between the big spheres. 
A rather similar phenomenon occurs when, in the pres- 
ence of say a hard planar wall, one has a single big hard 
sphere (solute) in a sea of small hard spheres (solvent). 

The concept of depletion force was first introduced by 
Asakura and Oosawai over fifty years ago in the context 
of colloid-polymer mixtures. Ever since, a great number 
of papers devoted to depletion interactions involving 
hard-sphere mixtures have appeared in the literature. 
The approaches have also been varied ranging from 
experiment,— density functional theory ^^ S i iiS t S t iSi com- 
puter simulations ; 11 ! 12 ' 13 ! 14 ' 15 ! 16 ' 17 virial expansions^ 
the Derjaguin approximation ; 19 ' 20 ' 21 ' 22 or the 
integral equation formulation of liquid state 

t ncor ., 7,23.24,25,26,27.28,29,30,31,32,33,34 

Very recently we have addressed the problem of deriv- 
ing the density profiles of multicomponent hard-sphere 
mixtures near a planar hard wall^ 5 . using an alterna- 
tive approach to the integral equation theory for the 
structural properties of the system i 36 ' 37 The main aim 
of this paper is to complement our former results with 
the study of the depletion potential for various systems 
involving hard-sphere mixtures. Specifically we will deal 
with the depletion interaction between two different (big) 



hard spheres immersed in a multicomponent mixture of 
(small) hard spheres. A limiting case of this situation 
is when the diameter of one of the big spheres becomes 
infinite and so this sphere is seen as a hard planar wall 
both by the other big sphere and by the multicompo- 
nent mixture of small spheres. It should be emphasized 
that our approach represents an improvement over the 
Percus-Yevick (PY) theory. In fact, it retains the major 
asset of the latter, namely it also yields analytical results 
in Laplace space. Further, our approach is technically 
simple too, so that the improvement over the PY theory 
is not achieved at the expense of adding difficulty to the 
theoretical development. 

The paper is organized as follows. In the next section 
we provide a relatively simple derivation of the Asakura- 
Oosawa depletion potential by looking at the exact re- 
sults for the radial distribution functions of a multicom- 
ponent hard-sphere mixture to first order in density. This 
is followed in Section IIIII by a summary of the results 
for the structural properties of a multicomponent hard- 
sphere mixture obtained by using the Rational Function 
Approximation (RFA) method. 36 ' 37 This section also in- 
cludes the limit where two of the species in the mixture 
(the solute) are present in vanishing concentration and 
the special cases where the solute particles have equal 
diameters or one is seen as a wall both by the other so- 
lute particle and by the solvent. Section [TV] presents the 
results both of the depletion potential and the depletion 
forces for some illustrative cases and compares them with 
available computer simulation data^ ' 12 ' 14 ' 17 We close the 
paper in Section |VJ with a discussion of the results and 
some concluding remarks. 



II. RADIAL DISTRIBUTION FUNCTIONS TO 
FIRST ORDER IN DENSITY: THE 
ASAKURA-OOSAWA RESULT 

In this section we present a simple derivation of the 
Asakura-Oosawa depletion potential that follows from 
exact results of the structural properties of multicom- 
ponent hard-sphere mixtures rather than the original 
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geometrical arguments. Let us consider an (N + 2)- 
component hard-sphere mixture in which the closest dis- 
tance between a sphere of species i and a sphere of species 
j is (Jij. Thus, <Ji = an can be considered as the "di- 
ameter" of a sphere of species i. The mixture can be 
either additive, i.e., aij = (er, + aj)/2, or nonadditive, 
i.e., <iij ^ {<Ji + <7j)/2. The number density of the mix- 
ture and the mole fraction of species i will be denoted 
by p and Xi, respectively. To first order in p, the radial 
distribution functions g^ (r) are exactly given by 

g iJ (r)=9l- ) (r)+pg^(r)+0(p 2 ), (2.1) 

where 

g$\r) = Q(r-a ij ), (2.2) 

N+2 
k=l 

x(<T ife + a kj - r) 2 0(r - a ij )Q(a ik + cr fei {2j*|) 

In Eqs. (|2.2p and (|2.4[) 6 (a;) is the Heaviside step func- 
tion. A derivation of Eq. (12. 4p can be found in Appendix 

H 

Now we assume that the mole fractions of species 
i = N + 1 and z = TV + 2 (here labeled as z = a and 
i = 6, respectively) vanish, so that the other species 
(i = 1,2, . . . , N) constitute the solvent. In that case, 
the depletion potential u a b(r) for the effective interac- 
tion between the solute spheres a and b is given by 

(3u ab {r) = - lng ah (r), (2.5) 

where (3 = 1/fceT, ks being the Boltzmann constant and 
T the absolute temperature. According to Eqs. (|2.1[) and 
(|2.2p . to first order in p (and for r > a ab ), one has 

/8ua6(r) = -pff£(r). (2.6) 
so that, taking Eqs. (|2.3p and (|2.4p into account, 

AT 

Pu ab (r) = -p—— V" (cr M + (Tfc - r) 2 [r 2 + 2(cr ai 
12r 

i=l 

+a bi )r - 3(a ai - a bi ) 2 ]Q(a a i + o- 6j : - r). 

(2.7) 

If both solute spheres are identical (a a i = cm, &ab = 
a a ), Eq. (|2.7p becomes 

N 

/3u aa (d) = -p—^x l (2cr a i -o a - d) 2 {d + cr a + 4a ai ) 

i=l 

xQ(2a ai -a a -d), (2.8) 



where we have defined the distance d = r — a a . If, fur- 
thermore, the ai interaction is additive, namely 2a a i — 
a a + Cj , then 

N 

/3u aa (d) = -pj^^Xi(ai -d) 2 (d + 3<7 a + 2ai)<d(<Ji -d). 
»=i 

(2.!)) 

This result coincides with the Asakura-Oosawa 
expressions^ 2 - 

We now go back to the case a ^ b, define z = r — a ab 
and assume that the ab and bi interactions are additive. 
In the limit en, — > oo the sphere b becomes a wall and Eq. 
([2~7| reduces to 



Pu wa (z) = - 




x(2z + <j a + Aa ai - ai) 

xe^ + ^-z). (2.10) 

Again, if furthermore the ai interaction is also additive, 

JV 

(3u wa (z) = -p- ^2 x i (°i ~ z f (2z+3a a +a l )Q (ai - z) . 

»=1 

(2-11) 

This result also coincides with the corresponding 
Asakura-Oosawa expressions^ 2 - 

Note that the validity of Eq. (|2.T[) to first order in p 
actually extends to any interaction among the solvent 
particles, including the so-called Asakura-Oosawa model 
(at = 0, a ai > a a /2, a bi > a b /2), i.e., only the solute- 
solvent (ai and bi) and solute-solute (ab) interactions 
need to be those of hard spheres. We must also point 
out that Eq. (12. 7p . while applying to first order in den- 
sity only, is quite general in the following sense: (i) the 
solvent may be in general polydisperse, (ii) the solute- 
solvent and solute-solute hard-sphere interactions are not 
necessarily additive, and (iii) the two solute spheres may 
have arbitrary sizes. We remark that, in general, the de- 
pletion potential u ab (r) is not a polynomial function of 
distance but a polynomial of degree four divided by the 
distance between the centers of spheres a and b. Only 
in the cases a = b [see Eq. (|2.8p ] and b — > wall [see Eq. 
(|2.10p ] does the potential become a polynomial (of degree 
three). 

The results of this section are exact but restricted to a 
low-density solvent. In particular, the Asakura-Oosawa 
potentials turn out to be purely attractive (with a range 
corresponding to the diameter of the solvent particles) 
and scale with the solvent density. Neither of these fea- 
tures remains as the density is increased. While it would 
be nice to have some measure of the error made in us- 
ing Eqs. (|2.9p or (|2.1ip in actual situations, there is un- 
fortunately no clear-cut way to estimate such an error. 
Instead, we note that, according to the qualitative dis- 
cussion performed in Ref. H3 . the entropic force grows 
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faster than the bulk density and becomes repulsive for 
distances on the order of half the diameter of the sol- 
vent particles. Hence, in order to account for these and 
other finite-density effects on the depletion interaction 
one must adopt a different strategy and resort to ap- 
proximations. In the next section we present our analyt- 
ical approach, which includes the PY approximation as 
a particular case. 



III. THE RATIONAL FUNCTION 
APPROXIMATION METHOD 

In this section we start by recalling the main aspects 
of the RFA method for multicomponent hard-sphere 
mixtures 3 ^ and refer the interested reader to our recent 
review paper— and references therein for details. 

As in the preceding section, let us consider an (N + 
2)-component fluid of hard spheres of diameters Oi and 
mole fractions Xi (i — 1, . . . , N + 2). Now we restrict 
ourselves to the additive case, but otherwise the density 
p is arbitrary. The packing fraction of the mixture is 
V = ( 7r /6)p(c 3 ), where 



N+2 



(3.1) 



denotes the nth moment of the size distribution. Ac- 
cording to the RFA ) 36 i 37 the Laplace transform Gij (s) of 
rgij{r) is given by 



Gij(s) 



2ns 2 



[L^.B- 1 («)].., (3.2) 



where Uij = (a, +(jj)/2 and L(s) and B(s) are (iV + 2) x 
{N + 2) matrices given by 



L ij (s) = L^+L^s + I^'^, 
B i: j(s) = (1 + as)Sij - Aij(s), 



(2) 2 



(3.3) 



(3.4) 



Aij(s) 



pXi 



(2) 

+ip {(7is)a i L\ j ' 



In Eq. 



i-xy 



(3.5) 



(3.6) 



By construction, Eq. (|3.2[) complies with the require- 
ment lim^oo se a ' ijS Gij(s) = finite. Further, the coef- 
ficients of s° and s in the power series expansion of 
s 2 Gij(s) must be 1 and 0, respectively. This allows us to 
express lJ ) and l_W in terms of lJ 2 ) and a: 



\m=0 



N+2 



4°/ = A + X'oj + 2X'a -XpJ2 x k <J kL 



(2) 
kj ' 



(3.7) 



N+2 



(3.8) 



fe=i 



where 



2tt 
1^' 



A' 



6irr) (a 2 ) 



(3.9) 



In principle, L' 2 ' and a can be chosen arbitrarily with- 
out violating any basic physical condition. In particular, 
the choice iff = a = gives the PY solution ! 38 ' 39 Since 
we want to go beyond this approximation, we will de- 
termine the coefficients lJ 2 ) and a by taking prescribed 
values for gij{oij) and the associated thermodynami- 
cally consistent (reduced) isothermal compressibility \- 
Hence, in our case, 



(2) 



(3.10) 



and a is found to be the smallest real root of an algebraic 
equation. 

Here we will take for gij(oij) the accurate extended 
Carnahan-Starling-Kolafa (eCSK3) approximatio n 35 ' 40 

( x 1 3r? (a 2 ) q-joj V 2 (5-2 V + 2 V 2 ) 

9ij{(Tii) = —r, + ^Z-fW)^, 



(o 3 ) *i. 



12(1-77)3 
rf{l + tj) ( (a 2 ) aiaj 



6(1 -n) 2 \(a 3 ) a. 



(3.H) 



which is thermodynamically consistent with the (re- 
duced) isothermal compressibility \ derived from 
Boublfk's equation of stated namely 



f 



677 (a) {a 2 } 



(I-77) 2 (1 - Ty)^ (a 3 ) • 
27 - 877 - 8rj 2 + 4?7 3 (ct 2 ) 3 



3(1 -v) 4 



(a 3 ) 



3\2 1 



(3.12) 



In the case where one of the species (say j = N + 2) 
becomes a wall (i.e., xn+2 — > 0, <jn+2 —* 00), Eq. (|3.11|) 
reduces to 



gwi(&wi) 



1 



1 - v (i-v) < ct3 > 



377 (a 2 ) 77 2 (5 - 2t7 + 277 2 

Oi + 



3(1 -t,) 3 



(a 3 )' 



4r/ 2 (l + r/) ({a 2 ) 



3(l-77) 2 V(^ 3 ) 



(3.13) 



A. Infinite dilution of species a and b 



k=l 



Now we assume that the mole fractions of species i — 
N + 1 and i = N + 2 (labeled again as i = a and i = 
6, respectively) vanish, i.e., x a — * 0, Xb — ► 0. In that 
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case, those species do not contribute to the total packing 
fraction or to other average values: 



N 



n < 3. 



(3.14) 



We assume that this is the case, even if the diameters of 
the spheres of species a and b are infinitely larger than 
those of the solvent species. 

The limits x a — > 0, x b —> imply that the last two rows 
of the (N + 2) x (AT + 2) matrix A defined by Eq. (|3~5]) 
vanish, so that the (N + 2) x (N + 2) matrix B defined 
by Eq. (|3.4|) has the following block structure: 



#21 



Bin 
B 2 n 



B 



Nl 



D 



NN 





V 



-Alb \ 
~A 2b 



-A 



Ni, 



-A 



Nb 



1 + as 
1 + as 



(3.15) 



Analogously, 







C\ a 




(B- 1 )^ • 


• (B _1 )2iV 


C 2a 


Gib 


(B- 1 )^ • 


■ (B~ )jViV 


CNa 


CNb 








1 

{l+txs) 














1 

(1+as) 



where 

Cia 



N 



E( B ~ r 



i k 



1,- 



(3.16) 
,N, (3.17) 



fc=i 



with a similar expression for C{ b . 

Insertion of Eq. (|3~T6]) into Eq. jSJJ gives G aa (s), 
Gab(s), Gbb(s), Gai(s) and Gbi(s) for i = 1, . . . , N, and 
Gij(s) for i, j = 1, . . . , N. The latter quantities refer to 
the A^-component mixture solvent and, as expected, are 
not affected by the presence of the solute particles a and 
b. The solute-solvent correlation functions G a i(s) are 



G ai (s) 



2irs 2 



EMb- 1 ) 

3=1 



1. 



, N. 



These quantities have been considered elsewhere^ in the 
wall limit a a — > 00. Here we want to focus on the 
solute-solute correlations G a b(s) in the presence of the 
Af-component bath. The result can be written as 



G a b(s) 



1 



-L a b(s) 



+ 



as [ 2irs 2 
e i(^-»>G ai (s)A ib (s) 



N 

i=i 



(3.19) 



where Eqs. (|3 . 1 7|) and (|3.18[) have been used. Equation 
p,19p is the main result of this section and readily allows 
us to get both the depletion potential u a b(r) and the 
depletion force F a b(r). They are given by 



I3u ab {r) = 



ln ^ 1 [G ab ( S )] 



(3.20) 



f3F ab {r) 



-0 



du ab {r) _ g' ab {r) 



dr g a b(r) 
C-^sGgbjs) - <Jabe~^ s gab(<J+b)} 

C-^Gabis)} 



1 

r 

(3.21) 



where £ -1 denotes the inverse Laplace transform oper- 
ator. In Eqs. (|3.20p and (|3.2ip it is understood that 
T > <y a b since both the potential and the force are of 
course singular in the region < r < <7 ao . We recall that 
the PY results are recovered by setting a = 0. 

When the A^-component mixture solvent becomes a 
pure fluid (i.e. <Tj = o\ for i = 1,...,AT), one has 



Li 



in and Ai 



XiAi for i,j = 1, 



, N, and 



Aia = XiA a and Aib = XiAb for i = 1, . . . , N, where 



A m {s) 



(p 2 (<7is)afL < f 7 l + ip 1 (a 1 s)a(L'{ 



2r(l) 



-(2) 
lm ' 



+ tpo(<Tls)(T 1 L [ { 

In that case, Eqs. (|3.18[) and (|3 . 1 9|) become 

L al {s) 



m = l,a,b. (3.22) 



G al (s) 



(3.23) 



Gab(s) 



2-ks 2 I + as - Ai_(s) 
L ab {s) [1 + as - A^s)} + L al (s)A b (s 



2-ks 2 (1 + as) [1 + as - Ai(s)] 

(3.24) 

In what follows we will consider the particular cases in 
which the sizes of the two solute spheres are the same or 
when one of the solute spheres has an infinite size so that 
it is seen as a hard planar wall both by the other solute 
sphere and by the solvent species. 



1. Case o a 



<7h 



Let us suppose now that the two solute particles are 
identical and, for the sake of simplicity, that the solvent 
is monodisperse. In that case, Eq. (|3.24|) reduces to 



G aa (s) 



e- a " s L aa (s) [1 + as - A^s)} + L al (s)A a (s) 



2-ks 2 (1 + as) [1 + as - A x (s)} 

(3.25) 

The second virial coefficient B 2 and the "stickiness" 
parameter r -1 associated with the depletion potential 
u a a{f) are evaluated in Appendix [Bj It is shown there 
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that the depletion potential predicted by the RFA in the 
colloidal limit a a j a\ — > oo is narrower and much deeper 
than that predicted by the PY approximation. In fact, 
the combination of depth and width represented by the 
stickiness parameter t _1 is divergent in the RFA and 
finite in the PY. 



2. Case (Tbjoa — > oo 

In the limit en, — > oo the solute particle b is felt as a 
wall by both a solvent particle and by the solute particle 
a. Before taking the limit <Jb — > oo, let us introduce the 
shifted radial distribution function 



L iw = nag W i{a W i), 



(3.35) 



(3.26) 
(3.27) 

(3.28) 



is the Laplace transform of j a b( z ) an d F' ab (s) = 
dr a f,(s)/ds. In the wall limit o\, — > oo, Eq. (|3.27[) yields 



r wa (s) = lim — e ^ s G ab (s) 



7a6(z) = 9ab{z + (Tab), Z > 0. 

In Laplace space, 

G ab (s) = e-°"-»" K6r Q 6(s) - r^(a)] , 

where 

/■OO 

r a &(s) = / dze- sz lah {z) 



1 + as 



^i OT (s) + e a " s G al (s)A iw (s) 

(3.29) 



2tts 



i=l 



where in the last step we have made use of Eq. (|3.19|) and 
have defined 

T ^ \ - r Lab ^ -J r \ - v Aib & 
L aw {s) = lim , A iw (s) = lim 



(Th— too 06 



IT6— >0O (7 fc 



(3.30) 

From Eqs. (JO}, l|3"3]) . (j3TTj) . (|3T5|) . and (j3~T(Jj) we get 

T r \ _T (0) , T (1) i T (2) 2 /o qix 



A iw (s) = pXj ^(^s)^ 3 ^ +(^i((7 l s)cr i 2 L-^ 



_l ( \ T (2) l 



L iw = A' - TraAp x j a J g w:j {a wj ), 

3=1 



iv 



(3.32) 



(3.33) 



i A A' A W 

L iw = - + ya-j - na-pcTi ^ Xj(rjg wj (a wj ), (3.34) 

3=1 



where i = a, 1, . . . , N in Eqs. ([333]) - ([335]) . The corre- 
sponding expressions for the depletion potential and force 
are 

(3u wa (z) = -ln7 lUQ (z) 

= -ln/T^r^s)], (3.36) 



/3F wa (z) = -0 



du wa (z) = 7L(f) 

7tua(2) 

C~ l [sT wa (s) - 7 wa (0)] 



(3.37) 



/:- 1 [r toa (s)] 

In the case of a monocomponent solvent (i.e., <7, = <7i, 
i > 1), one has Ai W (s) = XiA w (s) and Eq. (|3.29|) becomes 

1 I a „,(s)[l + as-Ai(s)]+L a i(s)A u) (s) 
" ,alSj ~™ 2 (l + as)[l + as- Ai(s)] 

(3.38) 

Of course, the same result can be obtained from Eq. 



IV. RESULTS 

In this section we illustrate the results that one ob- 
tains using our approach by considering some represen- 
tative cases. For simplicity, we will restrict ourselves to a 
monocomponent solvent so that <jj = a\ (i < N). With- 
out loss of generality we will measure distances in units 
of (7i and so the important parameters will be the sol- 
vent packing fraction r\ and the size ratios T, = cr b /a a and 
R = a a jo\. In Figs. [THE] we present the curves obtained 
using both the PY theory and the RFA approach as well 
as the corresponding simulation data.— i 12 i 14 i 17 

As can be seen from the figures, the RFA results cer- 
tainly represent an improvement over the PY theory in 
all cases for both the depletion force and the depletion 
potential, yielding in particular much better values for 
the well depth in the depletion potential. Our analysis 
will begin with the cases where both solute particles have 
the same size (£ = 1), namely those in Figs. HHHl The 
RFA results are clearly superior to the PY ones in the 
region a a < r < a a + \c\. For larger distances, however, 
the RFA and PY predictions are hardly distinguishable. 
When R = 2 (sec Fig. [I]) the oscillations of both the RFA 
and the PY curves are slightly dephased with respect to 
the simulation data. A similar behavior is exhibited by 
the density functional theory shown in Fig. 2a of Ref. 
[bH . Figure Q] also shows that the depletion force is a 
much more stringent quantity than the depletion poten- 
tial. In particular, the PY theory predicts a local mini- 
mum of F aa (associated with an inflection point of u aa ) 
at r ~ cr a + 0.2(7i, while both the RFA and the simulation 
data present a monotonic increase of F aa in the region 
a a < r < a a + a\. If R = 5 (Figs. [2] and [3]), the per- 
formance of the RFA with respect to the depletion force 
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FIG. 1: (Color online) Depletion force and depletion potential 
between two identical (big) hard spheres embedded into a 
solvent bath of (small) hard spheres as functions of distance. 
In this case, E = 1, R = 2, and rj = 0.3. Solid line: RFA 
approach; dashed line: PY result; circles: simulation data 
from Ref . [M 



becomes poorer as illustrated in Fig. [2] for ?/ = 0.116, but 
the theory is still able to capture even quantitatively all 
the features of the depletion potential for rj = 0.1 and 
i] = 0.2. For rj = 0.3, paradoxically in contrast to the 
case R = 2 of Fig. [T] it follows correctly the trend of 
the oscillations but otherwise overestimates the barrier 
height. Although not shown here and most likely related 
with the previous deficiency, also for rj = 0.3 the force 
starts to present features that seem not to occur in the 
simulations. We will come back to this point later on in 
connection with the hard planar wall limit £ — > oo. 

Now we turn to a more stringent situation, namely 
the case where the depletion effect takes place between 
a hard planar wall and a solute sphere, i.e., £ — > oo. 
In this instance, as shown in Figs. [4] and OH the agree- 
ment between the RFA results and the simulation data 
is also reasonably good. Particularly rewarding is the 
fact that, at least for r\ ~ 0.1, one gets a good perfor- 
mance even when R — 10 (Fig. [5]). Analogously to the 
cases with S — 1, the RFA strongly improves over the 
PY results for distances z < |<7i from the wall but both 
theories practically coincide for larger distances. Also, 
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FIG. 2: (Color online) Depletion force between two identical 
(big) hard spheres embedded into a solvent bath of (small) 
hard spheres as a function of distance. In this case, E = 1, 
R = 5, and rj = 0.116. Solid line: RFA approach; dashed line: 
PY result; circles: simulation data from Ref. 
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FIG. 3: (Color online) Depletion potential between two identi- 
cal (big) hard spheres embedded into a solvent bath of (small) 
hard spheres as a function of distance. In this case, E = 1, 
R — 5, and results are displayed for three values of rj. Solid 
lines: RFA approach; dashed lines: PY results; squares: sim- 
ulation data for rj = 0.1 from Refs. Q3 (open symbols) and 
(filled symbols); circles: simulation data for rj = 0.2 from 
Refs. 03 (open symbols) andU (filled symbols); triangles: sim- 
ulation data for rj = 0.3 from Ref. 0. 



the PY theory predicts a spurious local minimum of the 
depletion force near z — \(J\. It should be pointed out 
that in the planar wall limit one also starts to get a pe- 
culiar behavior, not shown in the figures, for relatively 
low densities (r\ — 0.2 if R = 10). This behavior, also 
shared by the PY theory, has similar features to the ones 
mentioned in connection with the poorer performance for 
some systems having E = 1, that is, the appearance of 
spurious local minima in the depletion forces and of in- 
flection points in the depletion potentials. While in the 
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z/a 

FIG. 4: (Color online) Depletion force between a hard planar 
wall and a (big) sphere in a background fluid of (small) hard 
spheres. In this case, E — > oo, R — 5, and results are displayed 
for two values of 77. Solid line: RFA approach; dashed line: 
PY result; circles: simulation data for r\ = 0.1 from Ref. H2I ; 
squares: simulation data for r\ = 0.2 from Ref. [T2I . 
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FIG. 5: (Color online) Depletion force between a hard planar 
wall and a (big) sphere in a background fluid of (small) hard 
spheres. In this case, E — > 00, R = 10, and r\ = 0.1. Solid line: 
RFA approach; dashed line: PY result; circles: simulation 
data from Ref. [3 



case £ = 1 within the RFA approach they may have their 
origin on the decreasing reliability of the contact values 
of the radial distribution functions with increasing size 
disparity R, when X — > 00 the features have to do with 
the fact that in the hard planar wall limit the radial dis- 
tribution functions both for the PY theory and in the 
RFA may become negative around the first minimum, 35 
which is clearly unphysical. 



V. DISCUSSION 

In this paper we have derived the depletion force and 
potential between two (in principle different in size) large 
spheres whose interaction is mediated by the presence 
of a multicomponent hard-sphere mixture (the solvent) 
composed of smaller particles. This has been done by 
using the RFA approach and the end results turn out to 
be completely analytical in Laplace space. One may say 
that the RFA approach not only retains the analytical 
character of the PY theory and the good performance of 
this latter at long distances, but it has some further assets 
as well. Hence, apart from the elimination of the ther- 
modynamic consistency problem, it is also able to correct 
the main drawbacks in the PY formulation in connection 
with the present problem, namely the poor prediction 
of the short distance behavior and of the well depth in 
the depletion potential and the non-divergent character 
of the stickiness parameter in the colloidal limit. Since 
it seems natural that the most relevant part of the de- 
pletion potential be the one corresponding to short dis- 
tances, the improvement over the PY result in this par- 
ticular region may be considered as a success of the RFA 
approach. But the fact that after such an improvement 
one can also cater for the (correct) long and intermediate 
distance behaviors represents another nice feature of the 
approach. It should be clear that, while for simplicity 
in the illustrative examples we have considered that the 
solvent is a monocomponent fluid, our development is 
far more general allowing us in principle to examine the 
same problem but with the solvent being a polydisperse 
hard-sphere mixture^ As far as we know, no simulation 
data for such a system are available and so a compari- 
son in this instance is not possible yet. In general, our 
expectation that the RFA produces reasonably accurate 
results both for the depletion forces and the depletion 
potential for low and moderate densities, provided the 
solute-solvent size ratio R is not too big, is fulfilled. 

We have already pointed out one of the limitations 
of the RFA approach (also present in the PY theory), 
namely the fact that in extreme conditions it may lead 
to unphysical (negative) values for the distribution func- 
tion ^/ wa (z), which in turn yield spurious features in the 
depletion interaction u wa (z). One technical point must 
be mentioned at this stage. It concerns the choice of 
contact values for the radial distribution functions of the 
mixture and the isothermal susceptibility. While here we 
have considered the eCSK3 contact values [see Eq. (|3.1ip ] 
and an isothermal compressibility \ that is thermody- 
namically consistent [see Eq. (|3.12[) ]. the RFA approach 
does not forbid the possibility of other choices. For in- 
stance, for high values of R one could instead take the 
simulation results or the ad hoc proposal of Henderson 
and Chan^ for the contact value J wa (0)- However, we 
have checked that, when using the empirical 7,^(0), the 
region where 7 wa (^) takes negative values does not dis- 
appear, although those values become less negative. We 
would expect of course that the more accurate the con- 
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tact values of the radial distribution functions and the 
isothermal susceptibility we use as an input, the better 
the performance of our development. This, however, re- 
mains to be assessed. 

Finally, we want to point out that in this paper we 
have restricted ourselves to the infinite dilution limit of 
the two solute particles. This has allowed us to equate 
the depletion potential with the potential of mean force. 
If the concentration of the solute is increased, this ap- 
proximation will cease to be valid. An important asset 
of the RFA is that it also yields analytical expressions for 
the direct correlation functions and the bridge functions 
of the mixture. These expressions could in principle be 
used for finite concentrations of the solute, for instance 
following the formulation of the depletion potential made 
by Castaheda-Priego et alM This we plan to do in future 
work. 



Acknowledgments 

We want to thank J. G. Malherbe, W. Krauth, E. 
Allahyarov, and H. Lowen for kindly providing us with 
their simulation data. M. Lopez de Haro acknowledges 
the partial financial support of DGAPA-UNAM under 
project IN-110406. This work has been supported by 
the Ministerio de Education y Ciencia (Spain) through 
Grant No. FIS2007-60977 (partially financed by FEDER 
funds) and by the Junta de Extremadura through Grant 
No. GRU07046. 



APPENDIX A: DERIVATION OF THE EXACT 
LOW DENSITY BEHAVIOR OF g tj (r) 

In a general mixture, the cavity function correspond- 
ing to the pair ij is defined as yijir) = e 13 ^^^ gij(r), 
where 4 l ij{ r ) is the interaction potential. To first order 
in density, 

VlJ ( r ) = i + pyg)(r) + 0{ P 2 ), yg ] M = E x *v$k(r)> 

k 

(Al) 

I 



Equation (|2.4p follows from Eqs. (|A4[) and (|A8[) . where 
it is assumed that Oij < Oi k + a k j f° r a U se ts k}. 



where 

vSiW= / dr'f ik (r')f jk (\r-r'\), (A2) 

fij(r) = e^^M — i being the Mayer function. In the 
case of hard spheres, fij(r) = — 0(r — cry), so that 

V [ S k {r)=V{a tk ,a k] -r), (A3) 



9% k (r) = V(a ik , a k3 ; r) 9 (r - <j l3 ) , ( A4) 

where V(Ri, R2] r) denotes the intersection volume of 
two spheres of radii R\ and R2 whose centers are a dis- 
tance r < R1+R2 apart. If r < Ri — R2 (where, without 
loss of generality, we have assumed that R\ > R2) the 
small sphere is entirely contained inside the large one, 
so that the intersection volume is just the volume of the 
small sphere, i.e., V(R\,R2]r) = ^--Rf- On the other 
hand, if r > R\ — R2, the intersection volume is the sum 
of the volumes of two spherical caps of heights h\ and /12, 
respectively, i.e., V{R\, R2] r) = v(Ri;hi) + w(i?2! ^2), 
where we have denoted by v(R; h) the volume of a spher- 
ical cap of height h in a sphere of radius R. Its expression 
is 

v(R;h) = -h 2 (3R-h). (A5) 
3 

It remains to obtain hi and I12 in terms of Ri, R2, and 
r. A simple geometrical construction shows that 

h 1 +h 2 = Ri+R 2 -r, R\-{Ri-hi) 2 = Rl-{R 2 -h 2 ) 2 , 

(A6) 

whose solution is 



(A7) 



(A8) 



APPENDIX B: EFFECTIVE STICKINESS OF 
THE DEPLETION POTENTIAL 

The second virial coefficient associated with the deple- 
tion potential u aa (r) is 

Bo = -2tt f drr 2 le- 011 "^ - ll 



, {R x + R 2 - r){R 2 - Ri + r) {R x + R 2 - r){Ri - R 2 + r) 

ill = , flo — . 

2r ' 2r 

The final result is then 

'^-Rl, r<R 1 -R 2 , 
V(Ri,R 2 ;r) = { ^(Jii + R 2 - rf[r 2 + 2{R 1 + R 2 )r - 3(i?i - R 2 ) 2 }, Ri - R 2 <r<R 1 + R 2 , 
0, r > Ri + R 2 . 

I 
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From here one can define the "stickiness" parameter— 
r- 1 = 4 



6B 2 

3 



7T£r 

12 d 

<j% ds 



G aa ( S ) 



1 + a a s 



(B2) 



s=0 



Making use of Eq. (|3.25[) one gets an explicit expression 
of t _1 in terms of the packing fraction 77, the size ratio 
<7 a /<7i, the RFA parameter a, and the imposed contact 
values g a a{o~a), ffai(oai), and <7ii(cti). In the special case 
of the PY approximation (a — 0), the result is 

[12(1 + 27?) + 3(5 + 477)o-i/o- a 



2(1 + 2i7) 2 

+6(1 - r])(ai/va) 2 + (1 - vfi^/aaf] <B3) 
In this approximation, the stickiness parameter T- 1 is 



lower bounded by 



6r] 



'coll 



1 + 27? 



(B4) 



In fact, this lower bound is the value in the colloidal 
limit a a /o~i — > oo. Moreover, the PY contact values in 
the colloidal limit are 



g a a(cr a ) 

9al{o-al) 

9n(o-i) 



3r? 



ai 2(1- r;) 2 : 
l + 2r? 

(1-77)2' 

i + Vg 

(1-77)2- 



(B5) 



In contrast to Eq. (|B4[) . the stickiness parameter pre- 
dicted by the RFA (a ^ 0) in the colloidal limit becomes 



~Ji = 12(cTi/(T a )a*5aa(o' aa ) + 6?/ 



1 + 3a*[l - 2g Q i(g a i) + 2rign(ai)] - 6o;* 2 g a i(g a i) [1 - (1 - 77)g Q i(cr Q i) 
l + 277 + 6a*77[l-2(l-7 7 )g 11 (ai)] 

I 



(B6) 



where a* = ajo\. Equation (|B6[) implies that, unless 
5aa(ca) ~ (cq/cti), the stickiness parameter diverges in 
the colloidal limit, i.e. 



T con = H{o-i/o- a )a*g aa {cF a ). 



(B7) 



The behavior g Qa (cT a ) ~ (oa/fi) appears in the PY the- 
ory. However, other theories (like the SPT, the BGHLL, 
and the one proposed by us in Ref. [44]) assume that 
gaa(o-a) ~ {o a /ai) 2 , while g aa (<J a ) ~ (^a/ci) 3 in our 
recent proposal^ and In g aa (cr a ) ~ (o a /o~i) according to 
Henderson and Chan 42, A simple geometrical argument 
shows that the divergence of t _1 in the colloidal limit 
is not an artifact of the RFA. The parameter r _1 essen- 
tially measures the area (in units of er^) below the curve 
r 2 [g aa (r) — 1] between r = a a and r — > oo. For large 
a a ja\ the range of g aa {r) — 1 is expected to be of the 
order of o~±. Therefore, the area can be estimated as 

r coll ~ -s^[9aa(0-a) ~ 1^1 ~ — 9 aa {(T a ), (B8) 

in agreement with the leading term in Eq. (|B6|) . To refine 
that argument, let us define the range £ of r 2 [g aa (r) — 1] 
as 



r 



d_ 

Or 

9'aa(°a) 



In [r 2 g aa (r)] 



2 

On 



(B9) 



9aa{o- a ) 

Taking into account the definition of G aa {s) as the 



Laplace transform of rg aa {r), we have 



On 



lim s 

s — >oo 



se^ s G aa (s) 



O a gaa(o a ) 

In the PY approximation, the result is 

I-77 



Ccoll — °~1 



2(1 + 2t?) 



(BIO) 



(Bll) 



in the colloidal limit. On the other hand, the RFA yields 
in that limit 



coll 



0"! 



gal(o- a l) O a 37/ 



9a 



1 

a* 



(a a) 01 1 - V 

O a 3?7 



.(B12) 



ai 2(1 - 77)5 

(1(1 

If 9aa diverges more rapidly than a a /a\, we get 

Ccoii = a* o x = a, (B13) 

which is generally much shorter than the PY result. In 
general, we have 



'coll 



1 r> gaa{°~a) , 
l^M tcoll, 



(B14) 



where /1 is of the order of 1. In the PY approximation 
we get /i = 1/3(1 — 77), while /j, = 1 in the RFA. 
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